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Abstract 

We investigate the level of fine-tuning of neutralino dark matter below 200 GeV in the low- 
energy phenomenological minimal supersymmetric standard model taking into account the newest 
results from XENONIOO and the Large Hadron Collider as well as all other experimental bounds 
from collider physics and the cosmological abundance. We find that current and future direct 
Dark Matter searches significantly rule out a large area of the untuned parameter space, but 
solutions survive which do not increase the level of fine-tuning. As expected, the level of tuning 
tends to increase for lower cross-sections, but regions of resonant neutralino annihilation still allow 
for a band at light masses, where the fine-tuning stays small even below the current experimental 
limits for direct detection cross-sections. For positive values of the supersymmetric Higgs mass 
parameter n large potions of the allowed parameter space are excluded, but there still exist 
untuned solutions at higher neutralino masses which will essentially be ruled out if XENONlt 
does not observe a signal. For negative untuned solutions are not much constraint by current 
limits of direct searches and, if the neutralino mass was found outside the resonance regions, a 
negative /^-term would be favored from a fine-tuning perspective. Light stau annihilation plays 
an important role to fulfill the relic density condition in certain neutralino mass regions. Finally 
we discuss, in addition to the amount of tuning for certain regions in the neutralino mass-direct 
detection cross-section plane, the probability to obtain some parameter value if the allowed model 
parameter space is chosen to be scanned homogeneously (randomized). 
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1 Introduction 



Recently, the XENONIOO collaboration has released new results after analyzing 225 live days of data 
taking. Limits on the spin-independent elastic Dark Matter (DM)-nucleon cross-section, a^^, have 
been increased by a factor of roughly four with 2.0 x 10~^ pb as the minimal value of the upper 
limit on a^^ at a dark matter particle mass of 55 GeV [1] . This leads to futher tests for dark matter 
models. 

Furthermore, the ATLAS and CMS collaborations have presented their analysis of more than 
5 fb~^ of 7 TeV, also including 8 TeV, data and claimed close to 5 local sigma level the existence of 
a Higgs boson with a mass of approximately 125 GeV [21 [3]. This fact fits very well to the minimal 
supersymmetric standard model (MSSM) because its prediction of the lightest Higgs boson mass is, 
when the LEP limit is taken into account, between 115 - 135 GeV depending on the supersymmetric 
parameters, see e.g. '^UUUIUU^- 

The existence of DM is supported by various cosmological observations such as gravitational 
effects on visible matter in the infrared and gravitational lensing of background radiation. Its total 
abundance, that has important implications for the evolution of the Universe, has been precisely 
measured by the WMAP collaboration |TT] during the last decade. This requires that a different kind 
of matter beyond the Standard Model (SM) of particle physics must be postulated. One of the most 
popular and most intensive studied candidate is the so-called weakly interacting massive particle 
(WIMP) that may constitute most of the matter in the Universe. Cosmology provides therefore a 
good motivation for Supersymmetry (SUSY), since the MSSM possesses a natural WIMP candidate 
as the lightest supersymmetric particle (LSP) is stable due to i?-parity conservation [121 US] (for 
reviews see e.g. [TH [T5].). 

SUSY (for reviews we refer to |16 1 I17 1 [T5]) has moreover the ability to solve the famous hierarchy 
problem by introducing superpartners with opposite spin statistics to each SM particle such that the 
loop contributions from superpartners cancel exactly and the weak scale is stabilized. Since SUSY 
must be broken, however, these cancellations are not exact and the non-discovery of SUSY particles 
pushes the breaking scale further up. This separation of the weak scale and of the SUSY breaking 
scale raises the question how easily this stability can be maintained. We apply therefore in this 
paper a measure of naturalness \19\ [20j , which was used for electroweak symmetry breaking, to the 
dark matter sector and study the level of fine-tuning. 

There is a series of studies on dark matter in the framework of simplified variants of the MSSM, 
the so-called constrained MSSM (CMSSM), which possess universal supersymmetry breaking mass 
parameters at the grand unification scale (for example |2H [22t I23j). Due to the existence of the 
grand unification condition on the gaugino masses in these models, there exists a LEP limit on 
the lightest neutralino mass: they must be heavier than 46 GeV. According to reference |23j the 
lightest neutralino mass must be larger than about 200 GeV at 95 % C.L. after taking into account 
all relevant experimental constraints. 

Instead of this restricted class of models non-universal gaugino models within the framework 
of the phenomenological MSSM (pMSSM) (see for example [SI [251 |26l [271 128] ) have gained much 
attention. In the pMSSM low-energy input parameters are used with no high-energy relations 
between them. These models were used to explain the possible annual modulation signals of 
DAMA/LIBRA [29] and CoGeNT p] {e.g. [3ll[32l[33l[Ml[35l[36l[37l[38l [39]), as well as the 
excess of nuclear recoil events reported by CRESST |40j . This would be interpreted in terms of dark 
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matter with a mass between roughly 10 GeV and 30 GeV and spin-independent cross-section of order 
10-^ - 10"^^ pb. However, it was shown that Ught neutrahno dark matter scenarios consistent with 
DAMA/LIBRA, CoGeNT and CRESST within the pMSSM are disfavored by LHC constraints [41j. 
In addition there are discussions about the vahdity and natural consistency of these signals [l2] and 
XENONIOO [HllSlIll] (see also [i5]) as well as CDMS 06], since these experiments have excluded 
these "would be" dark matter signals anyway. 

We assume therefore that dark matter has so far not been detected and ask how natural or fine- 
tuned the left-over parameter space is. Specifically we study in detail the not so well investigated 
neutralino mass range less than 200 GeV, taking into account all collider, cosmological and fiavor 
constraints including the recent results of LHC Higgs researches as well as fiavor studies. Over this 
complete mass region we find valid scenarios that may have escaped every experiment so far. We 
will especially show that it is possible to fulfill the muon anomalous magnetic moment condition for 
positive gaugino masses and a negative supersymmetric Higgs mass parameter, the /i-term. 

This article is organized as follows: in the next section, we define the fine-tuning measures and 
fix our notation of the neutralino sector. Then, in section [3] the method of our numerical analysis 
and the SUSY parameter space is discussed. Our results will be presented in section |4] including 
discussions about the annihilation mechanisms for neutralinos, the mapping of the level of fine-tuning 
into the direct detection cross-section plane, the direct detection cross-section and its dependence on 
the sign of the //-term, how the muon anomalous magnetic moment can be obtained correctly with 
a negative //-term, and lastly about functional fine-tuning and the probability distribution. Finally, 
we conclude in section [H 



2 Definition of fine-tuning 



SUSY needs to be consistent with the electroweak sector of the SM and has to reproduce the correct 
Higgs and Z-boson masses when electroweak symmetry is spontaneously broken. The complete 
scalar potential reads as: 



y = (|//|2 + m%Ji\H'j' + \H+\') + (|/.|2 + mlj{\H',\' + \H-\') + [b{H+H- - H^^H^,) + c.c^ 
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where /i is the SUSY respecting Higgs mass parameter from the superpotential, and '^u^ mass 
terms of the two complex Higgs doublets Hu and from the soft SUSY breaking part of the 
Lagrangian, h the bilinear Higgs coupling and g and g' are the ?7(1) and SU{2) gauge couplings, 
respectively. Minimizing this potential gives the well-known relation for the Z-boson mass at tree- 
level: 
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where tan f3 is the ratio of the vacuum expectation values (VEVs) of the two Higgs doublet fields. 

From equation ^ one can immediately see that SUSY mass parameters of the order of the weak 
scale are preferred to avoid tuning of the Z-mass already at tree-level. To quantify this, we use the 
fine-tuning measure defined as the sensitivity of the Z-mass 119\ [50] : 
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The parameters pi that determine the Z-mass on tree-level are the two soft Higgs mass parameters 
{rriHu ^'^d ruH^) and the bilinear coupling h. We take the total measure of fine-tuning arising from 
these parameters as a summation in quadrature: 



Afot = 



with the individual Apj's obtained in reference [U] 
4^2 



A6 



1 



1 + ^ ^ tan^ 2/3 



tan^ 2/3 



COS2/3+ ^cos2/3 



nir 



1 



cos 2/3 H sin^ ^ 



nir 











m| 








At 









1 



COS 2/3 
1 



+ 9 ^ tan^ 2/3 



1 + 



COS 2/3 



+ 



m 



^ 9 ^ tan^ 2/3 
ml, 



(4) 

(5) 
(6) 
(7) 
(8) 



Here niA denotes the pseudo-Higgs mass. 



The quantity Atot serves as an indicator how well a specific SUSY scenario avoids an unnaturally 
large separation of the electroweak and SUSY breaking scales. Small values are favored since they are 
less tuned and are therefore viewed to be more likely than those with high values of Atot- Note that 
the contributions given above do not depend on the stop mass and, hence, the fine-tuning will not 
show and exponential increase with the light Higgs mass. This is different to the fine-tuning measure 
for the pMSSM used in reference ^28j who evaluated Atot including terms which arise at leading log 
level (LL). The problem with these additional contributions is their dependence on the chosen cut-off 
scale, such that their specific value is scale dependent. This is to some extend the arbitrariness of 
fine-tuning definitions and due to this reason we will not include LL or higher contributions here. 
Our cut-off scale is the SUSY breaking scale (we concentrate on low-energy SUSY) and therefore 
LL terms of fine-tuning disappear. 

Apart from the sensitivity of the Z-mass, a tuning of the light Higgs mass (defined analogously 
to equation Q ) has been discussed in reference [38] . In their set-up it has been found to be of order 
100 already, but, again, a scale dependence cannot be avoided. Further discussions of fine-tuning in 
the (C)MSSM may be found in 03 EOl ED |52] . 

In addition to this "parameter" fine-tuning an "equation-tuning" can appear when cancellations 
between different terms are a consequence of model specific relations. This will be the case for 
the direct detection cross-section of neutralino dark matter, a^^. To quantify this, we evaluate the 
sensitivity of a^^ analogously to the Z-mass tuning and discuss its implications shortly in section 
For completeness: 

with Pi = tan (3, Mi, M2, m^}, see also 
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Ml is the bino and M2 the wino mass parameter. 

To fix our notation, we give some details about the neutralino x? which is a mixed state of the 
neutral gauginos (bino B and neutral wino W^) and the two neutral higgsinos (down-type and 
up-type #°): 



NuB + Ni2W^ + NisH^ + Ni^H: 



(10) 
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The coefficients Nij = 1,2,3,4) are tlie components of the mixing matrix that diagonahzes the 
neutralino mass matrix: 



I Ml -M^sin6'iy cos/3 sin 6*^/ sin /3 \ 

M2 M z cos 9w cos P — cos sin/3 

— iVfz sin6'iy cos/3 Af^ cos 6*1^ cos /3 — 

\ sin 6*^1/ sin /3 — cos sin/3 — / 



(11) 



3 Numerical analysis of the parameter space 

As already mentioned, we study in this work the MSSM defined at the electroweak scale, the so-called 
pMSSM, with the eleven free parameters as in reference 



tan /3 , Ml , , M3 , Ma , , mj^ , nij^ , rriq^ ^ , 777.53 ' «o , (12) 

where tan/3 is the ratio of the VEVs of the two Higgs doublet fields, Mi (i = 1,2,3) the three 
gauginos masses, the CP-odd Higgs mass and /i the Higgs-higgsino mass parameter. We chose 
different masses for left- (mj^) and right-handed sleptons (mj^) but no distinction in generations. 
Furthermore, we assume the left- and right-handed squark masses to be degenerate but different in 
generations (m^j ^ and rriq^ ) . 

The trilinear terms are parameterized by in the following way: 

At = aoYtrrig^ , Ab = aoYbrUg^ , Ar = aoYr^mjjnj^ . (13) 

This implies that we use non-zero trilinear couplings that are proportional to the third generation 
squark masses for At and A^, not as in reference [M], and the geometric mean of the slepton masses 
for the leptonic trilinear term. Since the SM Yukawa couplings of the first two generations are known 
to be very small, we can safely neglect A^, A^, A^, i.e. we set them to zero. 



In order to calculate the level of tuning, the parameters of equation (12) are randomly varied in 
the following ranges: 

Ml G [10, 200] GeV, M2 G [100, 2000] GeV, M3 G [100, 4000] GeV, 

niA G [90,4000] GeV, E [90,2000] GeV, ao G [-4.0,4.0] , 



"^91,2 



e [400,4000] GeV, niq^ £ [200,4000] GeV, tan/3 E [2,65] , 

e [100, 4000] GeV, mj^ £ [60, 4000] GeV. 



An input parameter x is linearly scanned as x = Xuim + A (xmax — a^min), where Xmin/max is the 
minimal/maximal value of x, respectively and A is a random number in the range [0, 1] generated 
using Mersenne Twister [55^ [56] . 

We emphasize that we have studied not only the case n > but also /i < which will turn out 
to be important. One might wonder if we can manage to satisfy the limit of the muon anomalous 



magnetic moment with a negative /x-term. This issue will be discussed later in section 4.4 



For our analysis we take into account different experimental data from cosmology, flavor and 
collider physics. The dark matter abundance constraints arise from the WMAP analysis [Tlj of the 
determination of the relic density. We set the neutralino relic density, Qh'^, within the 2 a range [57] 
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Quantity 




Reference(s) 


nh'^ 


[0.089,0.136] 


[m 




(121.0, 129.0) GeV 


[2l[3t[58l[59l[60l[6Tl[62l[63] 


Br(B s-y) 


[2.89,4.21] X 10"^ 


[6i] 


Br (5, ^ 


< 4.5 X 10-9 


[65] 


Bt{Bu ^ Ti?) 


0.52 < Rbtu < 2.61 


[66] 


Bi{K pLv) 


0.985 < Ri23 < 1.013 


[671 




[0.34,4.81] X 10-9 


[68^ 


xixi) 


< 3 MeV 


[69 




< 100 MeV 


i69j 


Ap 


< 0.002 


m 



Table 1: The experimental constraints. 



O/i^ G [0.089,0.136], where experimental and theoretical uncertainties are included. The influence 
of the direct detection search - XENONIOO (2012) - of DM will be explicitly shown in the plots 
representing our results. Results of ATLAS and CMS analyses for a standard model-like Higgs boson 
mass [2l[3l[58l[59l|60l[6ll[62l|63] provide the allowed mass hmit: 121.0 GeV < ruh < 129.0 GeV. 
Since the recent results of both collaborations are still preliminary, we take the range not too 
restrictive. A smaller mass range would not affect our results except that the numerical simulations 
would be more time-consuming. Also the light gluino mass and light squark masses of the first two 
generations are excluded by the analyses of the signature of missing transverse energy \70\ [71] , see 
also [721 ESj) so we allow their masses to be greater than 800 GeV and 1 TeV, respectively. 

We also take into account the constraint coming from pseudo-Higgs boson searches [74^ [75] 
that have excluded a significant fraction of the — tan /3 plane at small Ma and large val- 
ues of tan /3. Moreover, LEP constraints are included in our study: the invisible Z-decay width 
r(Z — )• XiXi) < 3 MeV [69], the pair production cross-section a{ee — )■ XiX2,3) < 100 MeV [76] . 
Ap < 0.002 [77] and the mass limits of supersymmetric particles [M] [77]. 

The experimental constraints including flavor and collider physics applied in our analyses are 
listed in table [T| where Rbtu is the ratio between the SUSY and SM prediction of the branching 
ratio Br(i3„ — )• tv'), R123 the leptonic kaon decay quantity and the muon anomalous magnetic 
moment. For the latter quantity we use the 3 a range [68j. 

The supersymmetric spectrum is obtained by SuSpect ^78j using the default SUSY breaking scale, 
the neutralino relic density il/i^, the spin-independent cross-section with protons a^^ and the anni- 
hilation channels were calculated by the micrOMEGAs code |79[ 180], while the low-energy observables 
(Br{B — )• 57), Br(i?s — )• p^p~), Rbtu, Rm) and have been determined by Superlso \81\ . 

Note that we do not calculate a x'^ as a probability measure in this work but see every scenario 
that passes the constraints as equally probable. Only the envelope of our final plots, that mildly 
depends on the parameters, has a meaning: The area outside of it can never be reached with the 
chosen parameter space. We take the viewpoint that tuned scenarios are less likely to form valid 
models of the pMSSM. 
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Figure 1: The relic density and the different (co) annihilation mechanisms presented for negative 
/i (see color coding and table [2] for more information about the final states). We have applied all 
collider and flavor constraints. The black lines show the upper and lower bounds that will be applied 
on the relic density; G [0.089,0.136]. Here qq stands for a light quark and U for a lepton pair. 
The plot is similar for positive fi. 



4 Results 

4.1 Obtaining the correct relic density 

First we discuss the different mechanisms that bring the relic density into the cosmological interesting 
region. The relic density is basically set by the thermally averaged annihilation cross-section, {aannv), 
using a freeze-out mechanism. In figure [T] we plot the relic density on a log scale versus the neutralino 
mass as obtained by our simulation with micrOMEGAs (for simplicity we plot here Qh"^ < O.2.), where 
the most important annihilation channels are presented by the indicated color code. In general a 
very efficient mechanism is needed to end up with a high enough (aannv) to obtain the correct relic 
abundance. 

Lepton final states (dark-blue points in figure [T]) are the most dominant annihilation channels 
and are present in our complete neutralino mass range. We want to point out that in past studies 
of the pMSSM the scenarios due to light stau annihilations at neutralino masses between 60 and 
80 GeV have been missed. To obtain these a distinction of left- and right-handed slepton masses 
is necessary as we will discuss later. The two resonant Z- and light Higgs-boson annihilations can 
easily be seen at around 40 and 60 GeV (red points), respectively. Above a mass of roughly 80 
GeV the neutralino may annihilate into two VF-bosons (light-blue) . There are further scenarios with 
dominant annihilations into light quarks for > 80 GeV hidden behind the lepton final states. 
The lower branch (dark-green) corresponds to chargino coannihilations, but for our chosen range of 
the lightest neutralino mass this mechanism is too efficient and never produces enough dark matter. 
These coannihilations are therefore unimportant for our further discussion, as we are also taking into 
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[GeV] mj [GcV] 

Figure 2: Dominant contribution to the neutralino annihilation in the — a^^ plane. The required 
relic density and all collider physics constraints have been applied except the anomalous magnetic 
moment of the muon. Here qq stands for a light quark and ii for a lepton pair. For better visibility 
we display the results in two figures. 



account the lower bound from the WMAP measurement. On the other hand, slepton coannihilations 
will play an important role at neutralino masses above 90 GeV (light-green). They are accompanied 
by a region of top final states at approximately 180 GeV (orange). 

When mapping the different models with their mechanisms into the — a^^ plane, we find a 
rather well ordered picture (all constraints except the muon anomalous magnetic moment have been 
applied) . We show an example for a negative /U-term in figure [2j 

The quark final states at neutralino masses between 40 and 60 GeV are due to the s-channel 
Z-boson and light Higgs boson, h, resonances. The chargino mediated t-channel and neutral Higgs 
mediated s-channel annihilations into VF-bosons fill a band at a^^ ~ 10~^ pb (light-blue). In this 
region (red points above 80 GeV) one also finds quark final states via heavy neutral Higgs, H, 
and CP-odd Higgs, A, exchanges. The most prominent mechanism, the annihilation into a pair of 
leptons, is homogeneously distributed and mainly mediated through t-channel light stau or s-channel 
Z, h, A, H exchange and their interference terms. Slepton coannihilations are situated in an area 





initial state 


final states 


chargino coannihilation 


Xl xt 


sc, ud, tb, ez/g, ^z/^, ZW 


slepton coannihilation 


x?n 


7r, ThjWvr, Zt 




Xl en 


76 




Xl I't 






Xl I'e 


We, Zve 



Table 2: List of the dominant annihilation channels as obtained by our simulation. 
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Figure 3: The level of fine-tuning for our models in the m^j — a plane. The colored points show 
the fine-tuning Atot for some specific parameter combination. Points with Atot < 10 are given the 
value 10 and Atot > 1000 the value 1000. The lines represent the exclusion limits of XENONIOO 
(2011) and XENONIOO (2012), and the prediction for XENONlt, respectively. 



below the W^-bosons final states. The chargino coannihilations are not present as argued abov^ 

The arrangement of the decay mechanisms holds generally and shows, aside from a different scale 
of the direct detection cross-section, no difference between positive and negative /i (for more details 



see section 4.3 ) 



Note that we have calculated the annihilation cross-sections of two DM into two photons in our 
models that satisfy all experimental constraints except the muon anomalous magnetic moment, and 
we found that the order of {a^-)(^^^.yv) is far below the recent claim of a gamma-ray line in the Fermi- 
LAT data |82j . Thus, we conclude that a dark matter particle mass of order 130 GeV and a partial 
annihilation cross-section into two photons of approximately 1.3 x 10^^'^cm^ s~^ is not compatible 
with our models. 



4.2 Fine-tuning and the spin-independent elastic WIMP nucleon cross-section 

We will discuss the spin-independent elastic WIMP nucleon cross-section as function of the neu- 
tralino mass with respect to the fine-tuning measure, equation Q. In this section, we show our 
results applying all experimental constraints discussed above except the muon anomalous magnetic 

^Some scenarios of chargino coannihilations have been present even after cutting the relic density, however, these 
were extremely rare (of order 0.02%) and we removed them from our study. 
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Figure 4: Plots motivating our upper bounds on the /i-term (left panel) and (right panel). The 
strong increase of the minimal amount of fine-tuning with fi and the negligible decrease of a^^ with 
rriA are visible. They are shown for a positive //-term but are similar for negative values. 



moment. To stress its impact we will consider this quantity separately in section 4.4 

First, positive values of the fi-teim will be presented: In figure [3] low fine-tuned regions can 
be found near the Z- and /i-resonances down to 10"^'' pb and for a^^ > 10~® pb at masses above 
80 GeV. A great part of the latter region is already excluded by the current XENON limit, so that 
low fine-tuned dark matter preferably appears for positive fi at masses between 20 and 60 GeV. 
Moving towards the exact Z-resonance a rise in Atot can be observed, which can be understood as 
follows: In order to obtain the correct relic density the neutralino-neutralino-Z-coupling has to be 
decreased to compensate for the resonant enhancement of {(Jannv)- This coupling is determined by 
the higgsino components of the neutralino [83] : 

{aannv)oc\Ci^^\^o,{Nf,-N!,f . (14) 



\ \ z 

Increasing fi decreases C^^ ^ (see equations (24) and (25)) and reproduces the correct relic abun- 



dance close to the exact resonance but at the same time increases the fine-tuning (compare the left 
panel of figure |4]). Note that at the Z-resonance /i-term is fixed due to the relic density condition and 
the fine-tuning stays roughly constant when moving towards lower direct detection cross-sections. 

Highly fine-tuned models appearing in figure |3 the otherwise low fine-tuned regions show a small 
tan/3 and a large niA- In that case Ab, equation (6]), becomes very large. 

The models lying out of the Z- and /i-resonances show a clear tendency towards higher fine- 
tuning for smaller a^^ (see also [53j). To decrease a^^, the higgsino component of the neutralino 
needs to be reduced which in turn is achieved through increasing the /i-term. Away from the Z- and 
/i-resonances a small neutralino-proton cross-section can therefore for positive n only be obtained 
at the cost of higher fine-tuning. To decrease further a^^ one could, in principle, increase /x above 2 
TeV and rriA above 4 TeV, but the effect on a^^ is small and the fine-tuning becomes unacceptably 
large (see figure [4| . 

Now, we will turn to negative values of the //-term (see figure [s]). The big difference is that 
the parameter region has just begun to be probed by direct detection experiments. Compared to a 
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Figure 5: The level of fine-tuning for a negative //-term in analogy to figure |3j 



positive fi-teim the neutralino-proton cross-section is shifted to smaller values, so that every mass- 
region still offers models that explain the dark matter riddle with low electroweak fine-tuning. An 
interesting region occurs for light dark matter near the Z- and /i-resonances where the fine-tuning 
stays small even if a^^ is decreased to tiny values. 

Across the complete neutralino mass region, scenarios with large fine-tuning show up at a^^ w 
10^^^ pb below the expected XENONlt "exclusion line". These come from light stau mediated 
annihilations whose stau mass is strongly influenced by the off-diagonal elements of the slepton 
mass matrix (see figures [2] and [6]) . The stau masses are given by: 



m~ +m~ - -Ml cos 2/3 + 2ml 

IL IR 2 



II 



mj + (-- + 2sin^6'vi/)Af|cos2/3 j + 4m2(A^ - //tan/3) 



(15) 



There are two ways to get light stau masses: (i) mj^ ^mj^ and {ii) rrij^ ^ ^T^' ^^^^ case we 



can approximate as: 



m~ ~ ?7i~ H — — |m^(j4T- — //tan/3)| ~ -M^ — |m^(j4T- — /itau/ 



(16) 



where we used — 5 + 2sin^^ty ~ and assumed in the last step mj^ ~ Mz- If a large //-term is 



present the ri mass can be suppressed (even below 100 GeV) when (Ar — ^ tan /3) — ^ ~ iQ^ GeV 
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Figure 6: Models with very light staus and their level of fine-tuning. We show the cases for a positive 
(left) and negative (right) /i-term separately. 



The same is true for the second case. As an approximation for the stau masses we find here: 

2 2 m^(^,-^tan/3)2 
mi:, m~ H ^ . (17) 

One can see that the stau mass is primarily determined by mj and for "^^ ^^lem'^'^'^ ~ 100^ GeV^ 

'l 

the mass is suppressed by the off-diagonal elements. This explains why the models with the lightest 
possible staus are severely fine-tuned. 

In figure [6] one can also see that in the neutralino mass region ~ [60, 80] GeV exclusively 
light staus are present. Here, light stau annihilation is the only possible mechanism to obtain the 
correct relic abundance which was not considered in past works on the pMSSM [211 [28] . For these 
scenarios we find that the second case, i.e. mr ^ rrir , is dominant to produce small m=, . 

Below neutralino masses of around 35 GeV we find for both signs of /i a separation of high and 
low fine-tuned models. Scenarios with a small /x-term, i.e. small fine-tuning, posses a significant 
higgsino component so that next to the slepton annihilation an annihilation via the Z-boson and 
their interference term is dominant for {aannv). For the points with m^^ < 100 GeV only stau 
mediated neutralino annihilation is important and we find that mj ^mj . A detailed study on light 
staus may be found in reference [54] . 

At the /i-resonance, ~ 60 GeV, one can find scenarios with high fine-tuning at a^^ ~ 10"^*^ pb 
above the XENONlt "exclusion limit" (see figure [5]). In contrast to the Z-boson resonance, the /i- 
term for /i-resonant neutralino annihilation is not fixed through the relic density condition, because 
the Higgs-fermion-fermion coupling is dependent on a SUSY parameter; in the case for annihilation 
into hh the coupling constant is proportional to cos/3. Models that have a small value of tan/3, 
i.e. a rather large cos/3 and, hence, a large Higgs-ft-ft-coupling, need a high ;^-term to decrease the 
neutralino-neutralino-Higgs coupling in order to keep the overall neutralino annihilation rate fixed 
at the value that reproduces the correct relic density. Despite their tiny higgsino component their 



direct detection cross-section, o"^^, is not minimal. See section 4.3 for further discussions on this. 
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Figure 7: The level of fine-tuning and the direct detection cross-section for negative (left) and 
positive (right) values of fj,. 

(These models can also be found in the vertical red stripe at 60 GeV in figure [6]) 

Note that the LHC phenomenology of the neutralino mass range under 70 GeV has been studied 
in reference |84j . They find that the invisible h decay width may exclude parameter regions of 
small n and Mi. Taking the most conservative upper bound of BRinv < 0.65 they present how this 
constraint cuts into the low fine tuned regions for a positive //-term. The implications on our study 
are not significant. Besides, for a negative /u-term the constraints are not competitive at all. 



4.3 The direct detection cross-section and the /i-term 

From figure [TJ one can immediately notice a clear difference for the spin-independent neutralino 
proton cross-section when the sign of fi is fiipped. A negative value allows for much lower cross- 
sections for a given absolute value of fj, than a positive sign. We want to point out again, that the 
red points showing up at small absolute values of fj, have a large rriA and a small tan f3. 

Let us discuss in detail why the sign of the /i-term plays an important role for the SI cross-section 
calculation. The dominating terms in the SI cross-section come from t-channel light and heavy Higgs 
boson exchanges. The formula can be found in references jl4| , 185 1 [5^: 



Fhh FhIh 

1 ' 2 



(18) 



where Gi? is Fermi constant and mred is the neutralino- nucleon reduced mass: 



rurcd = • (19) 

m~o + rriN 

The functions Fh^u and Ih,H are defined as follows 

Fh = {-Nil sin Ow + Ni2 cos ew) (A^is sin a + Nu cos a) , (20) 

Fh = (— A''ii sin6'iy -I- A''i2 cos^H/) (A''i3 cosa — A^i4 sina) , (21) 

Ih,H ^ 5^A;^^'^m,(iV|gg|iV) . (22) 
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Figure 8: The dependence of a as a function of the quantity Mi sin 2/3 and the wino component 
(^12) as obtained by our simulations. 



The angle a is the mixing of the mass eigenstates {h and H), and the coefficients kq'^ are given by 

^«-type = COS a/ sin p , k'^.type = - sin «/ cos /3 , 
ku-type = - sin a/ sin /3 , k^_^y^^ = - cos a/ cos (3 , 

for the up-type and down-type quarks, respectively. We neglect threshold corrections for simplicity 
in our discussion. 

To find a good approximation for the elements of the neutralino mixing matrix, we use the 
large SUSY scale approximation (Mj it |/u|)'^ ^ {i = 1,2) from reference |87j and find for the 
components: 

iV. = _Mlcos«.si„.„.^^^^^!«^ , ,23) 

. ^ Ml COS /3 -|- // sin /3 
iVis ^ -MzsmOw M2_l^2 ' (24) 



. ^ Ml sin /3 -|- u cos /3 

7Vi4 ~ MzsmOw \X 2 ^ • 25 

M-f — fj.'^ 



The unitary condition on the mixing angles yields 



iVii = Vl-^^i'2-^^i'3-^i'4 • (26) 

These approximations are in good agreement with the simulations. The only strong deviation 
occurs for scenarios with a dominant annihilation into VF-bosons near neutralino masses of around 
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80 GeV. For all other annihilation mechanisms they give sufficient precision for a reliable qualitative 
discussion. Evaluating the complete expression in the decoupling limit[^ yields the following formula: 



vr 



m 



red 



M| sin 



^/x COS 2/3 + ^ (Ml + ^ sin 2(3) 



1 2 



m 



H 



mi 



{NiismOw - Ni2Cosew) ■ (27) 



Note, that cos 2/3 is negative, i.e. cos 2/3 ~ —1 (tan/3 > 2), such that both contributions within the 
square brackets seem to have a different sign. However, looking at the ^-coefficients one can see that 
Ih and Ih preferably have opposite signs in the decoupling limit, such that both terms actually add 
up for a positive /U-term. 

On the contrary, if /i is negative, cancellations between both terms within the square brackets 
are possible and a^^ can be significantly smaller. However, this is only correct when Mi + /u sin 2/3 



is positive. Note that in this case the wino component, equation (23), is negative since M2 > Mi 



and \n\ > Ml for the great majority of our models (see figure [s]). The boundary at Mi + sin 2/3 ~ 
200 GeV occurs for the maximal values of Mi and tan/3 when, at the same time, is small. 

For Ml + fi sin 2/3 < a cancellation is no longer possible and the resulting a^^ is higher. Most of 
the scenarios with a positive wino component in figure [8] correspond to the earlier mentioned highly 
fine-tuned models at the Higgs-resonance. These showed large absolute values of the ^u-term and 
small tan/3. Then, Mi -|- /u sin 2/3 takes the smallest possible negative values and the cross-section is 
rather high {i.e. 10^^*^ pb). Also, those scenarios with light staus whose mass is suppressed by the 
off-diagonal mass matrix elements possess a positive wino component and hence map accordingly 
into the direct detection plane. 



An expression similar to equation (27) has been found in [88], however the wino component has 
not been considered there. Even though cancellations that occur between A^n and A^i2 may be 
neglected, since the bino component is in general much larger than the wino component, the sign of 
A'^12 is crucial for the behavior of a^^. 



4.4 The muon anomalous magnetic moment and the fi-term 

In figure [9] we applied the muon anomalous magnetic moment constraint {i.e. the deviation from the 
SM expectation) to our models with a positive //-term. In general it is not difficult to fulfill the 
condition. Compared to figure [s] the results essentially do not change. In reference [89j we find the 
three most important MSSM loop contributions to a^: 



f^Ti u ~ \ 3^ m^M2Mtan/3 , p \ 
af,{W-H,u^) = 4 , (28) 




a,{B,J.,-Jln) = ^^4^^^h^'^h (29) 



S7r2 Mf 

^5 u~ \ m^Mi^tan/3 ^ 



^I.e. we use sin^(Q: — /3) ~ 1, such that sin a ~ cos /3, cos a ~ — sin/3. 
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Figure 9: The level of fine-tuning after inclusion of the muon anomalous magnetic moment for 
positive values of the //-term. 



Here, the positive defined functions Fa and Ff, are given by: 

G3(x) - G^iy) 



Fb{x,y) 
Gsix) = 



X — y 
G4(x) - G4(y) 



1 



2{x - 1)3 

1 

2(x - 1)3 



X -y 

[(x-l)(2;-3) + 21nx] , 
[{x - l){x + 1) - 23;lnx] . 



It is commonly believed that cannot be fulfilled for a negative /i-term if Mi and M2 are posi- 



tive [90j. However, if equation (30), the bino-higgsino-right-handed smuon loop, dominates over the 



sum of the other contributions, equation (28), the wino-higgsino-muon sneutrino loop, and equa- 



tion (29), the bino-left-handed smuon-right-handed smuon loop, the total amount of the achieved 



positive pull of can be sufficient to correctly deviate from the SM prediction. This situation 

In this case the ri's are always light 



occurs in the limit of mj^/mj^ ^ 1, see left panel of figure 



10 



(< 400 GeV) due to equation (17). It is very important to note that is strongly dependent on 



the smuon parameters. Satisfying becomes easier once one abandons the slepton mass genera- 
tion universality since the relic density condition, which restricts the stau mass, would then become 
independent of a^. 



In our set-up we found that for large negative values of /i the contributions of equations ( 28 ) and 
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Figure 10: Only the green scenarios satisfy the condition. The left plot shows the condition 
as a function of M2/M1 and mj /mj . The latter ratio must be much greater than one in order 
to respect for a negative /i-term. The right panel indicates the behavior of when moving to 
smaller values of /i. 



(29) become important again due to the behavior of the loop functions Fa and Ff,. Because of the 



negative sign of // the pull of a„ tends to be toward negative values, out of the range that respects the 



experiments. We therefore find a lower limit on of about -1500 GeV, see right panel of figure 10 
This immediately favors low fine-tuned points as applying the condition automatically cuts out 
large negative values of /i. A great part of the models whose stau mass is strongly influenced by 
{At — fj,tanf3) drop out, and so do the highly fine-tuned scenarios from Higgs-resonant neutralino 
annihilation. Comparing to figure [2] one can see that a great part of the direct detection plane at 
low a^^ is filled by slepton annihilation. 

The muon anomalous magnetic moment condition has no effect on the low fine-tuned models 
and, as argued above, less fine-tuned scenarios are even preferred. A great part of the parameter 
space with Atot ^ 100 will be probed by future direct searches, but there are regions left that will 
not be tested. 

We emphasize that a negative sign of is by no mean in contradiction to the condition even 
though we did not distinguish between slepton generations. Therefore, this case should be paid more 
attention to in future studies. 



4.5 The functional tuning of a 



SI 



In section 4.3 we argued that the direct detection cross-section may become very small for a negative 
/i-term if cancellations between the light and heavy Higgs contributions are almost exact. This is 
only possible when the wino component is negative, or equivalently, for positive Mi -|- sin 2/3. In 
other words: This depends in a very sensitive way on the functional dependence of a^^ on this 
relation. 

To quantify these extremely accurate, tuned cancellations that might appear in cj^\ we therefore 
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Figure 11: The level of fine-tuning after inclusion of the muon magnetic moment for a negative 
//-term in analogy to figure [9) 



define a functional fine-tuning, Af, analogously to equation Q: 

aincr^l 



Afi^ 



(9 In Pi 



Af = 



E 



Pi=fj,,ta.n /3,-Mi ,M2,mA 



(31) 



The overall fine-tuning measure that includes both, the electroweak and the functional tuning, is 
then defined as: 



(32) 



In this way we deal with both tunings on equal footing. 

we show how the overall fine-tuning measure, ^fTj maps into the direct detection 



In figure 



12 



plane. Of course, this measure only makes sense for a negative /i-term since no cancellations are 



possible if /i is positive. It is visible comparing to figure 11 that for those scenarios in which the 
cancellations suppress a^^ below approximately 10~^^ pb the functional tuning measure becomes 
important and rises the overall tuning into an unacceptable range. For larger cross-sections the fine- 
tuning is still dominated by the sensitivity of the Z-mass and no change to the previously discussed 



results is observed. Comparing the left- and right-panel of figure 12 one can again observe the 
removal of the highly fine-tuned band containing models with light staus whose mass is strongly 
influenced by the off-diagonal stau mass matrix elements, when the constraint is applied. 

We see that even when the functional tuning measure is taken into account scenarios that possess 
a fine-tuning Uft lower than 100 are possible in the pMSSM avoiding all of our applied constraints. 
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Figure 12: The distribution of the fine-tuning measure ^ft, that includes the functional tuning 
of o"^^ and the electroweak fine-tuning, in the direct detection plane for a negative //-term. The 
constraint has been applied in the right panel. 



4.6 Probability distribution of our models 

Let us now consider the probability distribution of our simulated models in the direct detection plane 
which has so far been neglected in our discussions. We have shown several plots from global scans 
that covered the complete — a^^ parameter ranges as explained in section [sj In other words: 
We asked so far only how much tuning is required to reach a certain point in the — a^^ plane. 
It is obvious that the input parameter space is not covered homogeneously and we have avoided a 
statistical interpretation of the parameter mapping to the — a^^ plane. Note that the envelope 
makes sense and regions lying outside may never be reached. 

It is nevertheless instructive to look at a homogeneous (random) scanning of the input parameter 
space and to ask how likely certain regions in the — a^^ plane are reached. This is shown in 
figure [Tsj where the color coding shows in addition the previously discussed fine-tuning measure. 
It is immediately apparent that some regions are statistically more important. For example, we 
see that the regions from Z,h,H and A resonant annihilations and the light chargino mediated 
annihilation into VF-bosons have a stronger weight than light stau annihilations for positive fj, (see 



left panel of figure 13). Note, that the light chargino and heavy and CP-odd Higgs annihilations 
have now been ruled out by the XENONIOO (2012) update. For a negative //-term stau annihilations 
are also important, especially in the region G [65,80] GeV. Signals from direct searches should 
be expected in the untuned, densely populated regions of the parameter space, which are not yet 
excluded by data (especially XENONIOO). 



5 Discussion and conclusions 

In this paper we have analyzed the naturalness of neutralino dark matter in the non-universal 
gaugino model within the framework of the minimal supersymmetric extension of the Standard 
Model. We have taken into account all cosmological (upper and lower bound on the relic density), 
collider and flavor constraints including the results of XENONIOO (2012), LHC data on the Higgs 
mass, Br(i?s — )• and pseudo-Higgs searches. Hereby the soft supersymmetry breaking terms 
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Figure 13: The level of fine-tuning and the density of model points arising from a homogeneously 
(randomized) scan of the parameter space. The case of > is displayed in the left and < in 
the right panel. All experimental constraints are taken into account including the muon anomalous 
magnetic moment. 



are parameterized by 11 independent free parameters, that we have chosen such that the lightest 
supersymmetric particle is the lightest neutralino with a mass smaller than 200 GeV. We studied 
from the dark matter perspective how much fine-tuning is needed to reach a certain point in the 
— a^^ plane. Therefore we use a parameter fine-tuning measure (see equation Q) which was 
used before in order to study the naturalness of the Higgs - SUSY breaking scale separation. 

We first presented in figure 1 the contribution to the dark matter abundance Q.h? for different 
annihilation mechanisms as a function of m^. Demanding that neutralinos provide the right amount 
of dark matter we restrict the further scans to cases where nh^ £ [0.089,0.136]. 

We also investigated the dominant neutralino annihilation mechanisms (figure [T]) and have shown 
their arrangement in the m^— c^^-plane (figure [2]). In case of a signal in a direct detection experiment 
the most important annihilation mechanism can be deduced and we would know which channel at 
the LHC is promising for production of neutralino dark matter. To avoid the limits of the direct 
detection results of XENONIOO (2012), we showed that light stau annihilation of neutralinos in the 
early universe play a special role, not only in the mass range of light neutralinos, < 30 GeV, but also 
between ~ [60, 80] GeV, which has been missed so far in other studies. It is important to differently 
parameterize the soft SUSY masses of the left- [mj^) and right-handed {"mj^) sleptons. 

Note that we scanned the input parameter space in such a way that we obtain the fine-tuning 
for every point in the rriy^ — a^^ plane which is accessible. This implies that the density of points 
has no meaning, but that the envelope implies that certain areas cannot be reached by any input 
parameter. With this method we showed in figure [3] how the electroweak fine-tuning maps into 
the direct detection plane and found that a great part of untuned regions is already excluded by 
the current XENONIOO (2012) limit when the supersymmetric Higgs mass parameter, the /u-term, 
is positive. A general trend for higher fine-tuning for smaller a^^ is then visible. This trend can 
easily be understood, since the cross-section departs more and more from its natural value set by 
the generic scale. Thus, future direct detection experiments will push the amount of the electroweak 
fine-tuning further up. The only exception occurs for neutralino masses that allow for resonant 
annihilations, especially near the Z- and /i-boson resonances. 
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This last statement is also valid for a negative /i-term and the electroweak fine-tuning near the 
Z- and /i-resonance stays small independent of the value of a^^ is. Additionally, due to cancellations 
between contributions from light and heavy Higgs exchanges the direct detection cross-section gets 
shifted to smaller values, such that the XENONIOO (2012) exclusion limit is fulfilled easily. A 
negative //-term is therefore favored from a fine-tuning perspective (see figure [5]) . 

In our analytical study (see figure [S]) we have discussed the reason why the negative value of 
the /x-term allows for these cancellations in the spin-independent cross-section, and showed that the 
combination of input parameters {M\ + [x sin 2/3) and the sign of the wino component are responsible 
for decreasing a^^ to values that can be as low as 10^^*^ pb (see also figurejs]). Since these cancellations 
might be viewed as an instance of tuning, we reevaluated the fine-tuning by adding a measure of 
"equation-tuning" and found that scenarios with a^^ < 10~^^ pb are always unbearably tuned 



(figure 12). We find the scenario with the lowest fine-tuning, independent of the sign of at ~ 
84 GeV and cj^^ ^ 2.0 x 10"^ pb just below the new limit. 

Even for a negative /i-term we were able to get the correct positive pull for the muon anomalous 
magnetic moment, a^, to correctly deviate from the Standard Model (figure 10). This has been 



thought to be very difficult, but is possible due to the bino-higgsino-right-handed smuon loop which 
contributes to dominantly when both gaugino masses are positive (Mi > and M2 > 0) and 
mj^ ^ ^l^- If the latter condition is fulfilled staus are generally light (< 400 GeV) and help to 
respect the cosmological abundance of Dark Matter by light stau annihilation in the complete mass 
region of the neutralino. It should be stressed that there is an easy way to satisfy a^, namely, if 
additional parameters for the smuon masses, i.e. m^^^, are introduced. In this case we can avoid 
the connection between the relic density and the anomalous muon magnetic moment and fulfill 
both conditions without any doubt. Therefore, the case of a negative sign of the /i-term is equally 
important and should be investigated more carefully in future studies. 



Note that the density of points is meaningless except in figure [131 since we do not assign a 
probability measure, but determine only the amount of tuning required to reach a certain point 
in the — a^^ plane. The envelope implies, however, that these points cannot be reached. In 
this context it is interesting to note that the cross-section of the neutralino annihilating into two 
photons and into a pair of photon and Z-boson is loop suppressed and is therefore much smaller 
than the requirement from the claimed 130 GeV gamma-ray line in the Fermi-LAT data. Thus, our 
models can not explain this "evidence" . Besides, very light neutralino scenarios consistent with the 
DAMA/LIBRA, CoGeNT, CRESST experiments cannot be explained in the pMSSM especially due 
to the limits on Br(i?s — t- fi'^fj.^) and on the pseudo-Higgs mass-tan /3-plane. 



Finally, in section [461 we have also discussed in addition the probability distribution which arises 
from the the mapping of input parameters into the rriy — a^^ plane, (figure 13). We found that the 



Z- and /i-boson resonant areas become from both the statistical and from the fine-tuning point of 
views the preferred regions to detect neutralino dark matter if the /i-term is positive. For negative 
/i another statistically important region is formed by light stau annihilation, that has avoided direct 
searches so far. 
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